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Abstract 

In continuing the study of harmonic mapping from 2-dimensional Riemannian simplicial 
complexes in order to construct minimal surfaces with singularity, we obtain an a-priori reg- 
ularity result concerning the real analyticity of the free boundary curve. The free boundary 
is the singular set along which three disk- type minimal surfaces meet. Here the configuration 
of the singular minimal surface is obtained by a minimization of a weighted energy func- 
tional, in the spirit of J.Douglas' approach to the Plateau Problem. Using the free boundary 
regularity of the harmonic map, we construct a local uniformization of the singular surface 
as a parameterization of a neighborhood of a point on the free boundary by the singular 
tangent cone. In addition, applications of the local uniformization are discussed in relation 
to H.Lewy's real analytic extension of minimal surfaces. 

1 Introduction 



The classical Plateau Problem is the problem of finding a surface minimizing the area amongst 
all surfaces which are images of a map from a disk and spanning a given Jordan curve. We can 
formulate this problem more precisely as follows. Let A be the unit disk in R 2 . The area of a 
map a : A — > R n is 

A(a) = f 

The classical Plateau Problem. Given a Jordan curve T , let 

T = {a : A -> R n : a G IU 1,2 (A) n C°(A) and a\ dA : OA -> T is a homeomorphism] . 
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Find a* G T so that A(a*) < A(a) for all a E T . 

The classical Plateau Problem can be solved by finding a map a* which is conformal and harmonic. 
The study of classical minimal surfaces via harmonic maps is now well developed, starting with 
the solution of this problem by J.Douglas [D] in the 1930's. 

Our interest in this paper is a singular version of this problem. In |MYj . we investigated 
properties of 1-dimensional (2-dimensional resp.) subsets of a Euclidean space which minimize 
length (area resp.) amongst all the sets sharing a certain boundary set. Due to the shape of the 
boundary set, the minimizer is not a curve (surface resp.) but a union of curves (surfaces resp.). 
A variational formulation of this problem was presented in |MY] as a minimization problem for a 
certain energy functional of maps from a suitable polyhedral domain (or a finite complex). Con- 
tinuing the study of the two dimensional version of the problem, we refer to a singular surface as a 
union of the images of three maps, each from a disk joined along a curve 7 and bounding a graph 
T as pictured below. The boundary set Y in this picture is the graph consisting of three edges 
meeting at two vertices, and the curve 7 is the set where three disk-type surfaces are meeting. 
In the singular Plateau Problem, we look for a singular surface which minimizes area amongst all 
other singular surfaces bounding the given graph V . 




Figure 1: Singular Surface 

In order to capture the area-minimizer as an energy-minimizer in our singular setting where 
several (three) energy functionals are being minimized concurrently, we introduced in |MYj a 
weighted Z 2 -energy. As the name suggests, this functional is a quadratic combination of Dirichlet 
energies with respect to some weights. Its minimizer in the spaces of maps, weights, and conformal 
structures/glueing maps is shown to be the area minimizer. Furthermore, the minimizing map is 
harmonic and weakly conformal. The extension of the classical harmonic map theory to singular 
spaces was initiated by work of Gromov-Schoen |GS] and also Korevaar-Schoen |KS1] . [KS2j where 
singular target spaces are considered. The study of harmonic maps from a polyhedral space is 
considered in J. Chen [C], Eells-Fuglede (EE), Fuglede [FT] [F2] [F3] [F4] and Daskalopoulos-Mese 
[DMlj |DM2] [DM3] [DM4j . The use of the energy functional in singular minimal surface theory 
is a natural application of the singular harmonic map theory. 
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The best known singular surface minimizing the area is probably the (M, 0, <5)-minimal sets 
in R 3 (cf. |Alm] . JT]). J.Taylor [T] showed that in the interior of a (M, 0, <5)-minimal set, there 
are only two types of singularities. One is the Y-type, characterized by three minimal surfaces 
meeting at 120° degrees along a C 1,a curve (which we will call a free boundary). The other 
is the T-type, characterized by four of those free boundary curves meeting at a vertex, whose 
tangent cone is the cone over the regular tetrahedron and each free boundary is C 1,a up to the 
vertex point (0 < a < 1). In |MY] . we attempted to construct the singular configuration of three 
surfaces meeting along a free boundary as an image of a harmonic map from a simplicial complex. 
The key technical aspect in this approach is that there is an infinite dimensional moduli space V 
of conformal structures defined on the domain complex which is homeomorphic to three standard 
2-simplices glued along a common 1-simplex. Here the infinite dimensionality is caused by the 
different ways of glueing the three copies of the 2-simplices. The fact that the free boundary is 
real analytic in the (M, 0, #)-minimal setting suggests that the gluing maps one needs to consider 
for a minimizing element is smooth. 

The first result (Theorem [T]) in this paper is that when the conformal structure of a minimizing 
element of our variational problem is induced by a C 2 ' Q -smooth gluing map of 2-simplices, then the 
minimal singular surface has real analytic free boundary. It is an a-priori estimate, not unlike the 
result of Kinderlehrer-Nirenberg-Spruck |KNS] , in the sense that given a minimal configuration of 
a variational problem with a sufficient regularity, then the actual regularity is much better. Recall 
that |KNSj proves that the area minimizing singular surfaces of Y-type in R 3 which was shown 
to be C 1,a by [T] is in fact real analytic. The phenomena is accounted by the general principle of 
Lopatinski-Shapiro's, where an elliptic system of partial differential equations satisfies a coercive 
boundary condition. 

In order to set up our variational formulation according to this principle, we define a weighted 
^-functional of the Dirichlet energies. We show (Lemma [3]) that minimizing the ^-energy is equiv- 
alent to minimizing the / 2 -energy. The ^-energy minimizer satisfies the matching and balancing 
conditions (explained more precisely below, see (jlj) and (jSJ)) at the free boundary. This can be 
recognized as the coercive boundary condition in the sense of Lopatinski-Shapiro. 

Having established the real analyticity of the free boundary, we geometrically characterize the 
image of the minimizing map. Where the map is non-degenerate, namely away from the set S 
of isolated branch points, the image is locally either an embedded minimal surface or a union of 
three minimal surfaces meeting along a real analytic curve at 7r/3 angle (Theorem Ej). Note that 
this is consistent with the description of (M, 0, 5)-minimal surface by [Tj. Near a point on the 
free boundary parameterized by the non-degenerate map, the real analyticity of the three surface 
meeting along the real analytic free boundary is used to show that each of the induced conformal 
structures on the three minimal surfaces can be locally uniformized by a half disk, with a matching 
parameterization of the free boundary. This in turn says that there is a conformal (and harmonic) 
map locally defined on the unit ball (denoted Yq) about the vertex of the tangent cone at a point of 
the singular point. Thus we establish a local uniformization (Theorem [6]) of the minimal singular 
surface away from the branch points. 

Note that if the minimal surface is non-singular, this is nothing but the existence of an isother- 
mal coordinate system. The important difference in the two characterizations of minimal surfaces 
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via a conformal harmonic map is that in the classical Plateau Problem case, the parameterization 
is of the entire disk-type minimal surface, while in the singular case we only assert the existence 
of a neighborhood of a Y-type singular point conformally equivalent to Yq. This is because any 
two disk-type surfaces are conformally equivalent by the uniformization theorem which has no 
analog in the singular setting. From this viewpoint, the key observation contained in this paper 
and |MY] is that locally there is only one conformal structure Y to the singular surfaces while 
globally there are many, parameterized by the infinite dimensional space V. We further note that 
the conformal structure Y is conformally equivalent to the trivial element Xj^ (where the glueing 
via the identity map) of V . 

The existence of an isothermal coordinate system for the minimizer of our variational problem 
is closely related to the real analytic extension of each surface across the free boundary. The real 
analytic extension of a minimal surface across a real analytic boundary curve is a celebrated result 
by H.Lewy |Lej in the 50 's. The key to his proof is Schwarz' reflection principle of holomorphic 
function. Having the isothermal coordinate system around a Y-type singular point of the image 
of our minimizer, we construct the extension (Theorem [7]) of one surface via a weighted averaging 
of the harmonic maps parameterizing the three surfaces. We call this construction multi-sheeted 
reflection. The averaging works due to the important fact that one can regard these parame- 
terizations as a multi-valued harmonic junction defined on a half-disk and hence taking linear 
combinations preserves harmonicity. 

Furthermore when the ambient space is R 3 , the matching and the balancing conditions of 
the three minimal surfaces along a real analytic free boundary give an overdetermined boundary 
condition. Indeed one can construct the configuration from a single minimal surface with a real 
analytic boundary with the use of Bjoling's method. In this way, we provide an alternative 
expression (Theorem [8]) for the extended surface across the boundary curve. 

Acknowledgement. We would like to thank Joel Spruck for an explanation of the paper [KNSJ. 
The second author thanks the Japan-U.S. Mathematics Institute at Johns Hopkins University for 
the hospitality during his visits. 

2 Variational formulation 

In order to formulate the singular Plateau Problem precisely, we will need three copies of A which 
will be labeled Aj for i = 1, 2, 3. Let 

Bi = dAi n {y > 0} and A { = dA { n {y < 0}. 

We say $ = (0i,02 ; 03) is a gluing map if fa : dAi — > dAi is a homeomorphism satisfying the 
three point condition: 

fa(-l, 0) = (-1, 0), 0,(1, 0) = (1, 0) and 0,(0, 1) = (0, 1). 
We denote the set of all gluing maps by V . 
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We say $ G V is X-quasisymmetric (C k , C k,a resp.) if 0j is X-quasisymmetric (C fc (<9Aj), 
C k ' a (dAi) resp.) for all i = 1, 2, 3. For $ G P fl C 1 , we say $ is non-singular if 

for all p G 9Aj where # is the arc-length parameter of <9Aj. 

An extension of a gluing map $ G P is $' = (0' 1; 0' 3 ) where $ : A, — > Aj is a diffeomorphism 
so that 4>'i\ dAi = (pi for i = 1,2,3. We say an extension $' of $ is a Z)-quasiconformal (C k , C* ,M 
resp.) if 0^ is D-quasiconformal (C fc (Aj), C fc '^(Aj) resp.) for all z = 1,2,3. 

An element $6? defines a natural quotient space A$ of UAj as follows. First, we define an 
equivalent relation of UAi by setting s ~$ t if s G A^ t G Aj and 07* (<s) = ( t ) ] 1 {t)- We can then 
set 

X* = UA^ ~* . (1) 

The boundary dX<$> of Xq> is simply U5j. 

We note here a difference of treatment concerning X$ from our previous paper [MYJ, where 
X$ was defined using $ = (0 1; 2 , 03 ) where 0j's are quasiconformal diffeomorphisms of the disk. 
It turns out that X$ is completely determined by the parameterizations of <9Aj, as the different 
(quasiconformal) extensions inside Aj with the same quasisymmetric boundary map are regarded 
equivalent. 

Let T be a graph embedded in R n consisting of three arcs Tj (i = 1, 2, 3) sharing common 
end points g_ and q + . We will refer to such T as the fixed boundary of the Plateau Problem. 
We say \& = {ipi, 1^2,^3) is a uniform parameterization of r if ipi : — > R n is a constant speed 
parameterization of Tj so that g_ = ipi(— 1, 0) and q + = 0). 

Fix a uniform parameterization ^ of F . For each $ G P, we define JF,j,($) as the set of 
a = (oti, a 2 , "3) so that : A~ -> R n G WA 2 (Aj) fl C°(A~) satisfying 

a « 0«Ui = a j <Aj'U- (matching condition) 

and 

a i 0i| Bi = (boundary condition). 

The matching condition says that one can view a as a continuous map from Xq, and the boundary 
condition says that a restricted to dX$ parameterizes the boundary F as specified by \l/ and $. 
The matching condition implies that Oj(Aj) = aj(Aj) for i,j = 1,2,3. Let 7 = ai{A{). We will 
refer to 7 as the free boundary corresponding to a. Define the area of a as 

3 

A(a) =$^(oi). 
i=i 

We can now formulate the singular version of the classical Plateau Problem as follows: 
The singular Plateau Problem. Given a uniform parameterization ^ ofF, let 

t= U 
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Find a* G T so that A(a*) < A(a) for all a G T . 



A well-known method of obtaining the solution of a classical Plateau Problem is by the use 
of the energy functional. Motivated by this, we define a relevant notion of energy related to the 
singular Plateau Problem. First, we set 

C = {c = (ci, c 2 , c 3 ) G R 3 |ci + c 2 + c 3 = 1, q > for % — 1, 2, 3}. 

For a given c G C, we say that a G ^ r *($) is compatible with c if E(ai) — whenever Cj = 0. 
Otherwise, we say a is incompatible with c. We define the area of a as 

3 

A(a) = X)^(«<) 

and the ^-energy with weight c as 

/ 1 \ 1/2 
E (a) — \ ( E - (-^( a *)) 2 j ^ a ^ s compatible with c 

oo if a is in compatible with c 

where V] denotes the sum over i with q 7^ 0. 

c 

The Variational Problem. Find <£>* eV,c* G C and a* G satisfying the equality 

£ C »(V) = inf inf inf EJa) (2) 

The following theorem says that solving the variational problem above solves singular Plateau 
Problem. 



Theorem |MY] Suppose $* G V, c* G C and a* G satisfies (0|). T/ien a solves the singular 



Plateau Problem, i.e. 

A{a) < A(a) for all a G (J J r *($). 
Moreover, a* is a weakly conformal map and 

a- 



We have the following existence theorem. 
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Theorem (cf. |MYj ) Let V(K) be a subset of V which are K -quasisymmetric. For each K G 
[1, oo), there exists <& K G V(K), c K G C and a K G J r y($ K ) satisfying 

E cK (a K ) = inf inf inf EJa). (3) 

$ep(K)ceCaeT*(®) 

In particular, if 

inf inf inf EJa) = inf inf inf EJa), 
/or some sufficiently large K , then we can solve the singular Plateau Problem. 

Proof. Let $W g V(K),c (k) G C be a minimizing sequence in the sense that 

lim inf E r (k){a) = inf inf inf EJa). 

We recall that X-quasisymmetric maps have D-quasiconformal extensions for some D (cf. [?]). 
Let c&W be a .D-quasiconformal extention of By a slight modification of the proof of Theorem 
14 of [MY] , there exists a c^-energy minimizing map a^ G jF^,($( fc )), i.e. 

£? c C*o(a: (fc )) = inf E c(fc) (a). 

(We note that the modification of the proof of Theorem 14 only involves the observation that 
a i V 9 !^ has bounded energy independent of k since (p^ is D-quasiconformal.) Following 
the proof of Theorem 16 of |MYj . we see that (by taking a subsequence of k if necessary) 
$W _> $^ ;C ( fc ) — > c K ,a^ — > a A ' G .^(cl?^) and these limits satisfy (El), q.e.d. 

Our goal in the following sections is to investigate the minimizing map a* of ([2]) in the case 
we have an a priori C 2,a regularity along with a non-singular condition of the gluing map $*. 



3 Regularity of energy minimizing maps 

Fix $ G V and c G C. We say a G JF,j,($) is a c- weighted Z 2 -energy minimizing map if 

E c (a)= inf 

In this section, we prove higher regularity of a assuming a regularity condition on the gluing map 
$. In order to do so, we define another notion of energy which we call the ^-energy. 

Let k = (hi, k 2 , k 3 ) with ki > for i = 1, 2, 3. Given a G ($), the fc-weighted ^-energy of a 
map a is given as 

n 

E k (a) = Y / k l E(a l ). 
i=i 

We say a G JF,j,($) is a /c-weighted ^-energy minimizing map if 

E k (a) = inf E k (a') 
<*'e.F*($) 

We prove the following regularity result for an ^-energy minimizing map. 
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Theorem 1 For any k = {ki,k2,k^) with ki > and $ G V, let a G ^($) fre a k-weighted 
I 1 -energy minimizing map. If $ G C 2 ' M ; i/ien t/iere exists an extension $' 0/ $ so i/iai $' G C 2,z/ 
/or any < f < and a* o fa i is real analytic in A j U A, . In particular, oii o fa is a real analytic 
on Ai. 

Remark. Note here that such a gluing data $ in P belongs to V(K) for some K < 00. Hence 
the hypothesis of this theorem fits in the context of the existence theorem |MY] referred above. 

Proof. By |LT] (Theorem 4.3), there exists an extension fa { : A — > A of fa : S 1 — ► 5 1 
so that 0^ is a harmonic map from the Poincare disk (A,Oh yp ) to the Poincare disk, which is 
quasiconformal and fa { G C 2 ' u (Ai) for any < v < fj, . The harmonicity implies that fa i is real 
analytic in the interior Aj. 

Let A+ = {(x,y) G A, : y > 0} and U = {(x,y) G dAf : y = 0}. Let p : A+ -> A~ be 
the conformal map which fixes (1,0), (—1,0) and (0,1) and set <f>" = fa i o p. For simplicity of 
exposition, we will prove that u% = ai o fa[ is real analytic on Af U /j. This proves aij o fa = m o 
is real analytic since p~ l is real analytic o n A j U Ai. 

Let l g be the pull-back metric via fa[ : Af — > Aj of the standard Euclidean metric go defined 
on Aj. If we write l g = {gap) in terms of the Euclidean coordinates on Af C R 2 ( that is 
x 1 := x,x 2 := y,) then l g a/3 is real analytic in Af and C 1,u in Aj* - U J. Let 



H:={h = {hi, h 2 , h 3 )\hi : A+ — > R n G I^ 1 ' 2 }, 
and consider the functional £ :7i ^ R defined by setting 

<9/ij dh 



S(h) = J2h [ ■ ^-dx'dx 2 . 

2 = 1 I 



Set u = (mi,-u 2 ,-u 3 ) G Ti. Note that 0" is an isometry of (Af, % g) and (Aj,go), and thus 

£{u) = E k {a). 

This in turn implies that u is the minimizer for £(•) with respect to the elements h Eli, with the 
same boundary condition as u. By result of [M] and [DMlj . Ui is Lipschitz continuous in A^ U /j. 

We now claim C 1,13 regularity of Ui in Af U for /3 < |. For simplicity, we will use the 
following conventions. Let / = (fi, f 2 , ^3) be a set of functions /j : A+ U Ij — > R n . We set 
|£ = (ZL 9h d£\ and W = 9t2 8ta\ if 6£ £>& ex i s t f or i = 1,2, 3. We say / satisfies the 

ox \ ox ' ox ' ox J oy \ oy ' oy ' oy J ox 1 ox ' ' J J 

matching condition if 

/j(a;,0) = / J (a;,0)for i,j = 1,2,3 (4) 

and balancing condition if 

Now assume / satisfies the matching condition. We write / G C°° if /j G C°°(Af\Jli) for i = 1,2, 3. 
We write / G L 2 (L 2 oc , t^ 1 ' 2 , resp.) if /< G L 2 (A+) (L 2 oc (A+), ^ 6 2 (A+), W^(A+) resp.) 
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for i = 1,2,3. Furthermore, we set ||/||z,2 = X)f=i ll/«IU 2 with other Sobolev norms are defined 
analogously. We say / has compact support if h has compact support in A+ Ulj. We will also set 



Df = ^ max | 



dfi 



dx 



dfi 



and 



dy 



, max 



dx 



dh 



dy 



max 



dx 



dh 



dy 



f(x + h,y)- f(x,y) 



Claim. Define a functional O by setting 



M = E/ A+ v^§^ 

1=1 I 



(6) 



dx a dx 13 

for any rj G W 1,2 with compact support. If £ G C7°° with compact support, < ^ < 1, then 

\\CD{S h f)\\h < C (\\f\\wv + 0(<U(C%J))) (7) 
where C is a constant depending on | \ l g al3 \ \c^ a , k and (. 



Proof of Claim. In this proof, we use C to be a generic constant only dependent on 
iV^llc 1 .") k and (. We first replace rji by 8-hVi i R &■ Since 



we obtain 



E / V*(* + M) (**J£) ^ = -E / fcOfc + OOUi). 



1=1 



1=1 



Now replace ^ with ( 2 5hfi- If we let RHS and LHS be the right hand side and left hand side of 
©, we get 



< q|/|| wl>2 (\\5 h fD(\\ L * + \\((D5 h f)\\ L2 ) + 0(5_ h (( 2 6 h f)) 



(9) 



and 



LHS = f / fc, + K y) ( 5h ^L\ ( Sh W) C 2 + 2 / ^ ^(x + y, h) (s h ^\ c|W 



i=i 



' dx 13 



dx a I dx 13 



The last equality implies that 

/ dfA /, <)f. 



Y,Jkr 9 a ^x + h,y) k^j l^ h ^y 2 <LES + C\\(D6 h f)C\\ L ,\\DaS h f)\\ L ^ (10) 

Therefore, (jSJ), j9j) and (TX0|) combine to give 

\\C(DS h f)\\l 2 

< C(RES + C\\(D5 h f)(\\ L 2\\D((5 h f)\\ L2 ) 

< C(\\f\\ wl 4\\5 h fDC\\ L2 + \\C(D6 h f)\\ L2 ) 

+0(^(C 2 ^/)) + C||( j d^/)c||l 2 |I^C(^/)IIl 2 ) • (ii) 

Applying the inequality | \Shf\ \l 2 {q') < H-D/Hl 2 for any Q' CC X Id (cf. Lemma 7.23 |GT] ) and the 
Cauchy-Schwarz inequality 2ab < ea 2 + \b 2 several times in fTTTl) . we obtain 

\\CD{5 h f)\\h < C(\\f\\wi, 2 + 0(5„ h (( 2 5 h f))). 

q.e.d. (Claim) 

Denote the component functions of Uj as u™. Since the map u is ^-minimizing, each coordinate 
function u™ is also minimizing and the first variation formula gives 



3 

= " 

i=i 



for any r? G W ' with compact support. For any Qo CC Xj^, we can then apply Claim with C = 1 
in Qo,fi = uf and O equal to to conclude that ff^ G L 2 (Q ), Le. ^f- G W^' 2 



Now we choose rj G W 72 ' 2 map with compact support with compact support in Q CC X~ 



Id- 



Since G W 1,2 (f2), we can integrate the pointwise equality 



A I i g c*0 du T d Vi \ = (JL *g<A + i <& 9 ( ^It + i g c*i3 du ? 9 ( 9 lA 

dx \ dx a dx 13 J \dx J dx a dx 13 dx a \ dx ) dx 13 dx a dx 13 \dx J 

Note that the left hand side of the resulting integral equality is equal to since x t— > r)i(x,y) has 
compact support for all y. Furthermore, the third term on the right hand side of the resulting 
integral equality is equals by (fl2l since |^ G W 1,2 with compact support. Thus, we obtain 

9 i„ a 3 9 < 9 Vi ^ /",. i_afl 9 ( 9u ?\ 9 Vi 



hi K dx 9 d**dxf>-h] 1 9 



w., ,,, . ., i).r" \ dx 1 Ox 1 

i=i 1=1 \ / 

Let Q' CC Q and let C G C°° with < Q < 1, C = 1 in Q' and compact support in Q. We apply 

dx 



the claim with fa = -rf- and 0(rj) equal to the left hand side of the equality above to obtain 



du m , 2 / ' du m ' 2 



IIC^-^-)lll»<c^(||-^-||^ (n) + o(^ 
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Substituting rji with 5-h(( 2 ^h-^r) in 0(rj), we obtain 



am 



i=l 



^ ' ki dx 9 O.c" 



dx 13 



du T 



< C||u m || w i, 2 \\(5_ h 5 hir7r \\ L 2 + \\D 



< C\\u r 



\\(D(5 h 



du m 
dx 



\v + \\D- 



.du m 
dx 
du m 



\L 2 {Q) 



dx 



< C[-\\u 



m\\2 
\W 



-^ + -\\CD{5 h 



du m 
dx 



+ ~\\D 



L 2 (n) 
du r< 



~dV llL2 ^ ■ 



The above two inequalities combine to show 



\\(D(5 h - 



dx 



\h<c\ 



which in turn implies 



<2„,m 



d 2 u 



(dx 



>)<C\ 



du m 2 
\-fa\\w^(n) 

du m 2 

1-^-11^(0)- 



By the W 1>2 -trace theory, |^ G Wf ' 2 (A+ U /<) implies that G Lf oc (Ii). Therefore, for 

(a,0),(6,0)G4 



<9x 



(6,0) 



9a; 



(a,0) 



6 d 2 u m 
a dx 2 



x,0)dx < (b-a) 1 ^ ^jf 









dx 2 



for any g < 2 (and hence any p > 2.) This implies that u™ restricted to A is C 1,/3 for any (3 < 1/2. 
Thus, by the boundary regularity of harmonic maps, we have that G C 1,/3 (A+ U ij). 

We will now assume, for the moment, that u™ G C 2 (Af Uij) and prove that is rea l analytic 
m A+ U /i. The minimizing property of u™ says that is a smooth harmonic function in A+ 
(with respect to the metric l g) and satisfies the harmonic map equation 

d 







dx? 



^ da- 



Multiplying the above equality by ( = (Ci, C2, C3) £ C°° with compact support and integrating, we 
obtain 



i=i 



am 

■ a2 l -Cidx. 



dx a 



Since G C ,Q (Aj), we obtain the following linear system satisfied by 3n unknown variables u\ 



{El 



_d_ 

dx c 



dxt 



( B ?) 

(CD 2>'«r & „ 



UJ" - Mi = 

3 am 
in2 uu i 



J=0 fori = 1,2, 3, m = l,...,n 
fori = 2, 3, m—l,...,n 

= for m = 1, n. 



i=i 



11 



We wish to assign a weight the equations, unknowns and the boundary condition so that the system 
given above is elliptic and coercive as in [KNS]. We choose a weight s™ = to each equation (E™), 
t™ — 2 to each unknown variable -u™ - , r™ = —2 (i = 2, 3) to each boundary condition (13™) and 
r™ — — 1 to each boundary condition (Cj 71 ). The principle part of (E™), (resp. (B™), (C™)) is the 
part of exactly order 2 = s™ + t™ (resp. = r™ + for i = 2, 3, 1 = r™ + f™). Thus, we get the 
principle part of the linear system corresponding to these choices of weights as: 

d 2 u m 

{ET) ^8^ = ° for< = 1 ' 2 ' 3 ' "» = l.-.» 
(5™) M ™- M ™ = fori = 2, 3, m = l,...,n 

(CD E^y^ = 0, form = l,...,n. 
To show our system is elliptic, it is enough to show that any solution of (E™) of the form 

u m = c m e v^TK+OT) ? with o ^ ^ e R 2 

is trivial. From 
we obtain 

f V n crC 2 + { g 22 cTv 2 = o 

1 2 ^ 12 c-^ = 0. 

Since (£,77) 7^ 0, this implies c™ = 0. Hence (E™) is elliptic. 

To show our system is coercive, we need to show that any bounded solution (E™) of the form 

uT = e^ x v T(y), witho^eR 

satisfying is trivial. From (E™) we obtain, 

/)2 rre j to rl 2 tn m , 

in/3 _ iJ _ Ll i I i ll c 2, m , i 12 / , i 22 " SKj ^y^Tzg n 

and hence y?™ = c™ - e( a ^ +v/ ~ Tfc *"^ where z = a™ + >/— Tfr™ i s a ro °t °f the characteristic equation 

i„22^2 i o i„12 / Tt^ i„ll f 2 



y v + 2 y v-i£* - l g C = o. 

Since is bounded, a™ < for all i = 1,2,3 and m — 1, ...,n. The boundary condition 
implies 

and hence C :— C\ = ... — Cj. The other boundary condition (I?™) implies 

rrt 

fc=i 
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Thus, the condition that a™ < implies that C = 0. This shows u™ = 0. 

Now that we have shown our system to be elliptic and coercive with the chosen weights s™ 1 = 0, 
t™ — and r™ = —2 (k ^ 1) and r™ = —1, we can apply the elliptic regularity theorem of |KNS] 
(see also |ADN2| and [ME]). This implies that G C W (A+ U 

We are now left with proving that u™ G C 1,f3 implies u™ G C 2 . Indeed, the assertion follows 
from writing down the system {E™ 1 ) in divergence form then applying the Schauder estimate for a 
system in |ADN2j combined with a standard difference quotient argument (cf. [KNSj Lemma 5.1 
and the remark following it). Q.E.D. 

In order to prove higher regularity of an / 2 -energy minimizer, it is now sufficient to show that an 
Z 2 -minimizers is an / 1 -minimizer. We first need the following convexity and uniqueness statement 
for both the l l - and / 2 -energy minimizers. 

Lemma 2 Given two W 1 ' 2 maps a^^a^ : X$ — > R n in ^ r $($), define 
by setting af' := (1 — t)atj + ta^ . Then 

n . 

E k (a^) < (l-t)E k {a i0) )+tE k {a (1) )-t(l-t)Y,ki / \otf ] - af ] \ 2 dxdy. (13) 

i=i ^ 

If a continuous W 1 ' 2 map a : Xq, — > R ra zs k-weighted I 1 -energy minimizing in J r ^($), £/&en a is i/ie 
unique k-weighted I 1 -energy minimizing map. Similarly, if a is a c-weighted I 2 -energy minimizing 
in !F-%($>), then a is the unique c-weighted l 2 -energy minimizing map. 

Proof. The convexity of the Dirichlet energy (see for example, [H] or |SY] Chapter X (2.6ii)), 
E(af ] ) < (1 - t)E(ai 0) ) +tE(af ) ) - t(l - t) J \af ] - af ] \ 2 dxdy (14) 

immediately implies inequality (fl3|) . We note here that the cr^'s share the same Dirichlet boundary 
condition specified by \I/ and $. Thus, if and a'- 1 -' are both minimizing, then 

J |c4°' ) — a^\ 2 dxdy = 

for every i — 1, . . . ,n and this immediately implies the uniqueness of A;-weighted ^-energy mini- 
mizer. Furthermore, if af^ 1 ^ a\ , then Jensen's inequality applied to (|14|) implies 

E(a®) 2 <(l-t)E(a^) 2 + tE(aP) 2 
and this implies the uniqueness of the c-weighted Z 2 -minimizing map. q.e.d. 
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Lemma 3 If a is the k-weighted I 1 -energy minimizing map in then it is the c-weighted 

I 2 -energy minimizing map in with 

( E(a 1 ) E(a n ) \ 



\h T,j=i{E{a j )/kj) ' ' ' ' ' k n Y% = i(E(aj)/kj) J ' 

// a is the c-weighted l 2 -energy minimizing map in JF^,($) ; then a is the k-weighted l l -energy 
minimizing map in with 

k = (15) 

\ Xi x n J 

Here, x = (x\, ...,x n ) is the solution to the linear system (C — I)x = where C is a n x n matrix 
with all columns equal to c and I is the identity matrix. 

Proof. To see that a is c-harmonic, let a! be a competitor of a. We are assuming that 

n n 

£M?(ai)<£A^(50- (16) 

i=l i=l 

Use (fT6|) to show that 

n i 

E c (a) 2 = E-^(«0 2 
i=i Ci 

n i 

= V - E(a) 2 

^ E(oii)/(ki YJj=i{E{oij)/kj)) 1 l) 

n n i 

= ^hEi&^i-Eiaj) (17) 

i=i j=i k j 

n n i 

i=l i=l K 3 



it 



< Y-, 



(s(fii)A 



n i 

- ^E^/ik^m^/^f^^ 

n -, 
i=l Cl 



= E e (a) • E c (a'). 

Thus, we have shown E c (a) < E c (a') and this proves a is also c-weighted / 2 -energy minimizing. 
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To prove the last statement of the lemma, let at be the c-weighted Z 2 -energy minimizing, k as 
in (fl5|) and a' the unique /c-weighted ^-energy minimizing map. Since 



i=l 



the previous paragraph implies that a' is also the c-weighted / 2 -energy minimizer. The uniqueness 
of c-weighted Z 2 -energy minimizer implies that a = a' and hence a is also the fc-energy / 2 -energy 
minimizer. Q.E.D. 



Theorem 4 Let c G C and $ G T be so that $ G C 2,At and non-singular (i.e ^ 7^ 0). // 
a G JF^(<I>) zs a ^-energy minimizer with weight c, then there exists an extension G C 2,u of $ 
/or any < K /x so taat aj o fa g C w (Aj U Aj). In particular, a>i o fa is a real analytic function 
on Aj. Moreover, oti G C 1,7 (Aj U Aj) /or any < 7 < 1. 



Proof. The existence of $' G C 2 ' u so that that c^ o 
from Theorem [1] and Lemma [31 Furthermore, if we let tpi - 



G C w (Aj U Ai) follows immediately 
f>7 , then 



90 



(0) 



1 



Since ^-(^{(6)) 7^ this implies that ^ is Lipschitz. Moreover, 



difi 



09 



06 



1 



1 



1 



which immediately shows that -J^ is Lipschitz and hence 99^ G C 1 ' 1 . Thus, a 



atiO&ocfi G C 



1.1 



Boundary regularity theorem for harmonic maps then implies that ctj G C ' 7 (Aj U Aj) for any 
< 7 < 1. Q.E.D. 



4 Isothermal coordinates 

Let Af = {(x, y) G Aj : y > 0} and Jj = {(x, y) G <9Aj : y = 0} as in the proof of Theorem [IJ Let 
Y be the space which is the union of three closed half disks A+, Aj, A3 identified along the x-axis 
Ji, I 2 , 13 by the identity map. We denote the Ij's in Y by /. A map / = (/1, / 2 , / 3 ) : Y — > R n is 
said to be conformal if /j : Af — > R n is conformal (with respect to the standard Euclidean metric 
on A^) and satisfies the matching condition fi(x, 0) = fj(x, 0) for all — 1 < x < 1 and i, j = 1, 2, 3. 

Let Id = (idi, id 2 , ids) be the element of V so that idj : <9Aj — > <9Aj is the identity map. With 
this notation, let be as in ([!]) with $ = Id. We denote the Aj's in Xu by A. A metric on X^ 
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is g = { % g} so that % g is a metric on Aj. Let F = (Fi,F 2 ,F 3 ) : F - * (-^idj*?) be a map which 
is injective but not necessarily surjective so that Fi(Af) C A*. We say that F is an isothermal 
coordinate system at p G A if Fj(0,0) = p , Fj is a conformal map from A^ to (Aj, *g) and 
F(x, 0) = Fj(x, 0) for i,j = 1, 2, 3. 

Let re?,c'6C and a* G J r ($*) satisfy the equality © with $* G C 2 ' M and non-singular. 
Then a* is a weakly conformal harmonic map which in C 1,7 (Aj U A,j) by Theorem HI Thus, the 
differential da* can be continuously extended to be defined on A, U Ai. Let Si C Aj U A; be the 
set of isolated points where the da* is rank zero. We define the singular set of a* as 

S = {p G X$,* : p G Si for some 2}. 

and the regular set R to be X$*\S. 

Let p be a point in i? fl A Then near the point a* (p) , the image of the map is a collection of 
three embedded minimal surfaces meeting along a curve. This curve is real analytic since it is the 
image of real analytic map a* o (j>i\ A , by Theorem HI We can thus apply Corollary 12 of [MYj to 
see that there exists a neighborhood U of p so that three minimal surfaces a*(Ai n £7) meet along 
the free boundary e a * := a* (Ai fl U) at 120° angles. 

We summarize the preceding statements as follows. 

Theorem 5 Assume $* EV,c* G C and a* G satisfies the equality (TJ|). J/$* G C 2,a and 

non-singular, then for every regular point p G -R fl A C Xm, t/iere a neighborhood of the point 
a*(p), where the image of a* is an union of three embedded surfaces Si,S 2 and S 3 meeting along 
an embedded real analytic curve 7. The unit normal vectors 771 , 772, 773 of Hi, E 2 , £3 a/ong 7 satisfies 
the geometric balancing condition: 

X> = o ( 18 ) 
1=1 

ai every point 0/7. 

By reparameterizing Xj^, we have the following local uniformization of the singular minimal 
surfaces near a(p) for every regular point p G R fl A. 

Theorem 6 Assume &* E V,c* E C and a* G satisfies the equality (Tj^. J/$* G C 2,a and 

non-singular, then for every regular point p G RC\A C X^, i/iere exists conformal, harmonic map 
f = (fii J2, f'3) '■ Yq — ► R- n parameterizing the image of a in a neighborhood of a(p) 

Proof. By Theorem HJ there exists an extension $' of $* so that a* o <f> i is real analytic 
in Ai U Ai. Let % G be the the pull back of the Euclidean metric on R n under the map a* o 
and l g be the pull back of the Euclidean metric on Aj under the map <f>\. By Theorem HI l G is 
C w (AiUA,j) and *g is C 1,7 (Aj U Ai). As a c*-weighted / 2 -energy minimizer, a* is weakly conformal 
and hence a* o is weakly conformal with respect to l g. Therefore % G is conformally equivalent 
to l g, i.e. % G = \ l g for some nonnegative function A. With respect to the Euclidean coordinates of 
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Aj, denote the metric components of % g by l g a p, and % G by l G af 3. We now wish find an isothermal 
coordinate by solving the Beltrami equation 

w- = pW z 

where Beltrami coefficient p is given by 



p = 



5n + l g22 + U'gix l g22 - *<?i 2 ! Gn + 'G^ + 2/Gn 'G^ - *G 



i2 
r 12 



Note here that the Beltrami coefficient p is represented by the pull-back metric % G = (ajO0' i )*G ? o as 
well as by the pull-back metric l g = (</>9*<?o- As $ is a G 2 ' a -diffeomorphism of bounded dilatation, 
the Beltrami coefficient p has moduli strictly less than one and is C l,a . Furthermore, the metric 
components % G a p are given by 

Ga ?- { dx« ' dx? )R "- 
Thus, the components of l G are real analytic on A U Ai since the map o 0- is real analytic 



there. Note that the quantity y i G\\ i G22 — i G\ 2 is the pulled-back area form of the immersed 
surface o 0^(Aj U Ai) possibly with isolated branch points in R ra . Away from the branch points, 
namely in a neighborhood where the differential of the map o 0j is non-degenerate, the quantity 



G11 l C 2 2 — % G\ 2 is real analytic, as the term % G\\ ^22 — 1 G\ 2 is strictly positive. Hence /i is a 
ratio of two real analytic functions on the open set R where the differential of the map o i is 
of rank two as the latter function is strictly positive on the set. 

Recall here that there is the conformal reparameterization of the unit disk Aj by the half disk 
A^ via the Riemann mapping p : A+ — > Aj fixing (1, 0), (—1, 0) and (0, 1). Pull back the Beltrami 
coefficient /ij (which is covariant as p is conformal) and the regular set R defined on the disk, onto 
the half disks by the map p. By abuse of notation, we use pi and R for the pulled-back coefficient 
and the pulled-back regular set respectively. 

Let P = ai o 4>i(p) be a point on the free boundary with p = p(xo,0) G R H A. We set 
w(z,z) = u(x,y) + iv and p = r](x,y) + i((x,y) to rewrite the Beltrami equation defined on the 
half disk A + as the following system of equations with real analytic coefficients: 

( u v\ = ( C (I+V) Y 1 ( (l-V> x + (v x \ 
\v y ) \(l-rj) ( J \-{u x -(l+ri)v x J 

The inverse matrix on the right hand side exists because \p\ 2 = r/ 2 + ( 2 < 1. We also have the 
Cauchy initial data 

u(x, 0) = x — x and v(x, 0) = 

for (x, 0) G Ii near (xq, 0) G U. Therefore, we can apply the Cauchy-Kowalewski Theorem and 
obtain, in some neighborhood of the point p, a unique solution to the Beltrami equation. This 
solution Wi is a quasiconformal diffeomorphism from a neighborhood Ui C A+ U Jj of (x , 0) to 
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a neighborhood Vj C A+ U Ii of (0, 0). By construction, the pulled-back metric of the Euclidean 
metric go of A^~ under Wi is conformal to p*( l g). Thus the map w^ 1 provides a parameteriza- 
tion of the neighborhood of p _1 (p) = (xq, 0) in (A+ U Ii.p*^)) by a open set Vj in A+ U U. 
After scaling and then composing with p, we have constructed an isothermal coordinate system 
F = (Fx, F 2 , F 3 ) : Y -> (X u , g) of p e A. The map / = (fx, f 2 , f 3 ) with f i = ai o §\ o F satisfies 
the desired properties of the theorem, q.e.d. 



5 Multi-sheeted reflections via isothermal coordinates 

Let / = (fx, f 2 , i^)) : Yq — > R™ be as in Theorem [61 The equality 

/ i (x,0) = / i (x,0) (19) 

for i, j = 1, 2, 3 implies 

^(x,0) = ^-(x,0). (20) 
ox ox 



Using the conformality of fi, (118]) can be written as 



3 dfi 3 dfi 

o = E|r(M) = Etoo)- 

i=l I dy I »=1 I dx I 



This combined with with ([20]) implies 

n = ^ 

" dy 



= E^(^0). (21) 



If we let 

A(x, y) = -A(x, -y) + ^ E -2/) ( 22 ) 



3f =1 



then (USD and (J2U imply that 



f 1 (x,0) = f 1 (x,0) and ^(z, 0) = ^(x, 0). (23) 



We claim (T23l shows t/i : A — > R n defined by setting 

Ux(x,y) = 



fx(x,y) for y > 
/i(ar,j/) for y < 



is harmonic. Indeed, for any smooth £ : A — > R n with compact support, integration by parts 
gives 

- / V^-VUxdxdy = [ £Afxdxdy - [ £^(x,0)dx 
Ja+ Ja+ Ji oy 



and 

- / Vf • VU x dxdy = [ ^Af x dxdy + [ £^-(x, 0)dt 
J a- J a- Ji ay 

where A + = {(x,y) G A : y > 0}, A" = {(x,y) G A : y < 0} and / = {(x,y) G dA + : y = 0}. 
Summing up the above two equations and using the harmonicity of f\ and fi, we obtain 

- / V£ ■ VUxdxdy = 0. 

iA 

By Weyl's Lemma, is a C"^ harmonic map. Similarly, there exists extensions U~2, U3 of fi 
and /3. We call this construction of the real analytic extension C/j of fi the multi-sheeted reflection. 
By summarizing the argument above, we have 

Theorem 7 Lei £ 1; £ 2 , £3 and 7 as in Theorem^ The surface Sj can fre extended real analyt- 
ically across the curve 7. This extended surface is parametrized by the conformal, harmonic map 
Ui via the multi-sheeted reflection. 

We note that the extendability of the minimal surface Sj across a real analytic boundary curve 
7 follows from a celebrated result of H.Lewy [Lej . On the other hand, Theorem [7J gives a more 
precise picture of the extension. Indeed, the extension of the parameterization f\ of Ex is given in 
terms of a linear combination of odd reflections of fi, f 2 , fz as defined in (T221) . 

When three minimal surfaces are geometrically balanced along a curve in R 3 (i.e. the unit 
outer normal of the three surfaces sum to zero as in ffl8|) ). the entire configuration is completely 
determined by one of the three surfaces. This follows from the so-called Bjoling's problem resolved 
by H.Schwarz. We will explain below how to use this and arguments in the proof of Theorem [7] 
to give a construction of Lewy's extension. 

Theorem 8 A minimal surface in R 3 with a real analytic boundary can be extended across the 
boundary by a multi-sheeted reflection. 

Proof. We start with a surface Si with a real analytic boundary curve 7. Let rji be the unit 
outer normal to the surface £j along 7. Let n 2 and 773 be the two unit vector fields defined on 7, 
normal to 7, each making the angle of 7r/3 to T]i, Note here f]i + f]2 + V3 = 0- The solution by 
Schwarz of the Bjorlng's problem ([M] III §149) then provides locally defined, uniquely determined, 
minimal surfaces £2 and £3 along 7 so that 772 and 773 are unit outer normals to £2 and £3 along 
7 respectively. 

The next step is to construct real analytic local parameterizations of the singular surface 
(U 3 =1 £j)U7 by a continuous map from Y ; i.e. we construct maps from Afuli to £iU7 for i = 1,2,3 
so that their restrictions to U agree. For example, one can use the hodographic projection of |KNS] 
to devise maps satisfying the required condition. We briefly recall the hodographic projection. 
The surfaces £ 2 and £3 are locally graphs over the tangent plane T (? £ 1 . We suppose that £ 3 lies 
above the plane, £ 2 below. Let n 2 and LT3 be the orthogonal projection maps from £ 2 fl B$(q) and 
T, 3 nBs(q) to T g £i for sufficiently small 8. We can choose the coordinates so that q = (0, 0, 0) G R 3 , 
the tangent line to 7 at q is the xi-axis and T g £i is the x 1 x 2 -plane. Define U3 by the conditions 
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n2 1 (xi,x 2 ) = (xi, x 2 , u 2 (x)) G E 2 and (xx,X2) = (x\, x 2 , u 3 (x)) G E 3 . Near the origin, the 
map h(xi,x 2 ) := (xi,u 3 (x) — u 2 (x)) is of rank two. It sends n 2 (7) = 113(7) to the Xi-axis and its 
image is contained in the upper half plane. The map h o Yl 2 is the hodographic projection and its 
inverse map Il 2 1 o h^ 1 defined on a sufficiently small half disk centered at the origin defines the 
real analytic parameterization of E 2 around q, while IT3 1 o h~ l defines that of S3. Similarly, Ei 
can be parameterized using the tangent plane Tq(E 2 ). Denote these three maps parameterizing 
the neighborhood of q by H = (H u H 2 , H 3 ) : Y -> (Uf =1 Ej) U 7 with Hi : U k -> Ej. The real 
analyticity and the continuity of H follows from the construction. 

For each % = 1, 2, 3, the pulled-back metric H*Gq on Aj (recall that Go is the Euclidean metric 
on R 3 ) defines a conformal structure which in turn defines a Beltrami coefficient /ij. This Beltrami 
coefficient is real analytic due to the argument in the proof of Theorem [61 We again solve the 
Beltrami equation locally around a point p$ = (xq, 0) 6 A ; C via the Cauchy-Kowalewski 
theorem. The solutions to the three Beltrami equations give a parameterization F = (Fi, F 2 , F3) 
of a neighborhood of p in X id by (X Id , g ), an isothermal coordinate system. Define / = (fx, f 2 , / 3 ) 
by fi = 9i°Fi- Since /j : (A + , /*Gq) — > R 3 is an isometric minimal immersion, /j is harmonic with 
respect to the metric J*Gq (cf. [LaJ). As harmonicity is conformally invariant on two dimensional 
domains, /j is also harmonic with respect to the Euclidean metric on A+. Hence the resulting 
parameterization f\ : A^ U Jj — > (Ej U 7) C R 3 is harmonic, conformal and without branch point. 
Furthermore, we also have 



where c 



dfi 




dfi 


dx 




dy 



Now each fi : A^ — > R 3 can be extended across the real axis by f\ of 



In particular, we have a conformal parameterization of the extension of Ei as in Theorem[71 q.e.d. 
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